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Abstract — This  paper  intended  to  demonstrate  how  to 
regenerate  the  simple  forms  and  the  complicated  patterns 
composed  of  sine  curves  using  Fourier  series.  The  result 
shows  that  the  Fourier  series  is  a preferred  method  over 
other  methods  described  in  this  paper.  The  experimental 
result  of  the  regenerated  patterns  proves  the  efficiency 
and  the  reliability  of  the  proposed  method. 
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I.  INTRODUCTION 

Many  high  quality  and  high  trusted  printed  documents 
have  a set  of  forms  & models  constructed  upon  a family 
of  curves  and  specially  generated  using  limited  functions 
like  sinusoidal  formulas.  Some  of  these  curves  are 
boarders,  guilloches  and  a portrait  drawing  found  within 
the  banknote  document  and  stamp  document. 

Depending  on  the  characteristics  of  the  process  being 
studied,  a method  of  analysis  could  be  used  to  extract  the 
useful  data  picked  with  their  main  features  to  find  out  the 
hidden  relationships  and  implied  in  the  tabulated  original 
form. 

Different  mathematical  methods  may  be  used  for  estimate 
and  extract  the  suitable  featured  data  of  these  models. 
Fourier  series  method  can  be  used  among  others  like  least 
square  or  latent  Root  Regression  methods  to  do  the 
extraction  and  the  estimation  of  the  suitable  featured  data 
needed  for  the  reconstruction  & the  regeneration  of  the 
original  printed  models.  The  Fourier  series  and  Fourier 
transform  (FT)  is  without  doubt  an  important  tool  not 
only  in  the  field  of  modern  signal  - and  information 
processing,  but  also  in  other  engineering  sciences.  For 
instance,  Fourier  transform  contribute  greatly  to  the 
sampling  theory  [1],  time-frequency  analysis  [2],  and 
wavelet  theory  [3],  which  constitutes  the  foundations  of 
today’s  digital  world  [4,  5,  6].  One  of  the  important 
results  concerning  Fourier  series  and  FT  is  the  so  called 
Riemann-  Lebesgue’s  Lem  (RLL).  Its  version  for  the  FT 
asserts  that  the  FT  of  an  integrable  signal  / E L1  (R)  on 
the  real  line  has  a regular  behavior,  in  the  sense  that  it  is 
continuous,  and  vanishes  at  infinity,  i.e.  it  converges  to  0 
as  1 60 1 ->  oo  Analogously,  for  the  Fourier  series,  it  says 
that  the  Fourier  coefficients  of  a w.l.o.g.  27i-periodic 


functions,  which  are  each  integrable  on  w.l.o.g  [-71,  tt], 
decay  with  increasing  index  toward  0. 

II.  STATISTICAL  METHODS 

The  model  of  multiple  linear  regressions  explains  the 
dependent  variable  y;  by  a linear  function  of  k explanatory 
variables.  So  in  a sample  of  n observations  the  model  can 
be  written  as  [7] 

y = xp  + £ e (s)  = o 

E(e  £ ')  = <j 2 In 

(1) 

Where  x matrix  of  explanatory  variables  with  n x k,  y 
vector  of  dependent  variable  with  n x 1,  p vector  of 
parameter  with  k x 1 and  8 is  a vector  of  error  with  n x 
1. 

The  least  squares  estimation  (the  points  estimate)  of  P by 


minimizes  [7] 

s = (y-xb)  (y-xb) (2) 

So  ft  can  be  estimate  by 

P = (x'x)-1x'y  (3) 


The  matrix  of  x must  have  rank  k.  The  unbiased  estimator 
for  <J 2 is 


e'e 

(n-k) 


(4) 


Where  e is  a vector  or  error  with  nxl  e = ^y-Xpj  to 

measure  the  degree  of  linear  association  between  a vector 
y and  a set  of  k vectors  x we  use  the  coefficient  of 
determination  which  define  by  R2  with  equation  [7] 

R2=i-I>,2  /iL-y)2 (5) 

i=l  / i= 1 

So,  some  types  forms  for  simple  and  multiple  linear 
Regression  can  define  like  [7] 

1- 

yi=y?o  + Axi  + Ax2+-  + Axk  + £i 
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2- y.^Pqx/1  x/2 

3- Yi  =Po  +Pl  X1  +p2  X2  +-  + Pk  Xk  +ei 

4-  Yi  = — + Pi  x;+  £j 

a 

Some  basic  assumptions  for  model  (1)  are  not  satisfied  so 
we  can’t  use  equation  (3)  to  estimate  the  vector  p [7]. 

If  the  model  define  with 

y = x(3  + e E(e(  = O 

E(ee')  = Wct2 

(6) 

Where  W is  a sale  matrix  and  w"1  exist  the  estimates  for 
vector  p in  model  (6)  when  apply  the  method  of  least 
squares.  Define  by  [8] 

P = (x'  w"1  x)  1 x'  w_1y  (7) 


y = h(e)+s (12) 

Where  y is  n x 1 Vector  of  depending  Variable,  h(  0 ) 
is  a n x 1 Vector  consisting  a nonlinear  function  with 
independent  Variable  Xi  like 

h (x1 , 0) , h (x2 , 0) , , h (xn  , O)  , by  linearizing 

the  model  in  (12) 

y = h(e0)+^^(e-e0)+s  <i3> 

Then  the  approximate  relations  to  estimate  the 
parameter  depending  on  (GLS)  method  like 

0-00  + (z'z)_1  z's  (I4) 

Where 

7 an(e.) 

00' 


This  is  a generalized  least  squares  (GLS)  estimator  of 
vector  p [8]. 

In  nonlinear  regression  model  which  is  nonlinear  in  the 
unknown  coefficient  vector  p there  are  many  types  of 
nonlinear  model  like  [8] 

1-  The  univariate  nonlinear  regression  model  with 
function 

Yi  =h(xj  ,e)+£j  i = l, , n (8) 

Where  h is  a nonlinear  function. 


2-  The  M- Variate  nonlinear  regression  model  with 
function 


yij=hj(xi,0)  + fij  i = l,2, , n 

j = 1,2, , m 


(9) 


So,  there  is  M of  these  regression  equations  explaining 
n dependent  Variable. 


3-  A system  of  nonlinear  simultaneous  equations. 
The  model  define  by 


®j(yi>xj>0i)=eij  j = l>2, ,m 


(10) 


In  this  model  the  equation  j the  i th  endogenous  Variable 
is  an  explicit  function  of  other  endogenous  Variable  [9]. 


To  maximize  a function  L(  0 ) which  has  a 
continuous  first  and  second  derivatives  with  respect  to  the 
Vector  0 , there  is  many  Numerical  methods  of 
maximization  like  [10]. 

1-  The  gradient  algorithm  with  equation 

e,=e»+k„  u 1 (» . ) (is* 

Where 

f (0O)  the  gradient  evaluated  at  0o 
ko  a scalar 

Ho  matrix  of  information 


2-  The  method  of  steepest  ascent 

for  equation  (15)  Ho  = I (unit  matrix) 
ko  may  be  expensive  to  compute  since  it  requires  the 

0f  (0O) 

matrix  — rr: — of  the  second  partials  of  L with 


5 0' 

respect  to  0 

3-  The  Newton-Raphson  method  with 
n-i 

and  when  k = 1 the  equation 


H0  = 


01=0O 


0f(0j 
00' 
is  the 

3f(e0) 


which  gives  the  maximum  of  L(0)  because. 
1-1 


50' 


f(®o) 


.(16) 


For  model  (8)  the  function  Ln  to  be  minimized  define  as 

L(0)=^Z(yi-h(xi’0))2 00 

2n  tt 

For  simply  we  define  the  model  (4)  like 


4-  The  Gauss-Newton  method  is  an  approximation  to  the 
Newton-Raphson  method  in  that  when  the  matrix  of 
second  derivatives  of  L is  computed. 


www.ijaems.com 


Page  I 676 


International  Journal  of  Advanced  Engineering,  Management  and  Science  (IJAEMS) 
h ttps://dx.  doi.  org/1 0. 24001/iiaems.  3.6.10 


[Vol-3,  Issue-6,  Jun-  2017] 
ISSN:  2454-1311 


5-  The  method  of  scoring  the  Newton-Raphson  method 
for  maximizing  the  average  log  likelihood 

1 n 

— T^Log  L(yt  , 0)  uses  as  Ho 

n tt 


6-  The  method  of  quadratic  hill-climbing,  in  equation  (15) 
- Ho  can  be  replace  with 


gfjgo) 

ae' 


Where  a is  a scalar  is  chosen  to  maximize  L(  0 ),  the 
method  requires  computing  the  characteristic  roots 


of  the  matrix 


df(e.) 

5 0' 


III.  APPLICATION  PART 

Fourier  method  is  used  to  determine  the  coefficients  of  a 
Fourier  series  which  approximate  periodic  time  series 
data  which  can  be  defined  as  [11]: 

i=l,2, ,n 

Xy 

j=l,2, ,r 

Where  r is  the  number  of  Replication  and  n is  the 
number  of  observation  for  any  Replicate  so  Xij  is 


Replicated  values  for  ith  observation  in  jth  Replication  the 
formula  for  Fourier  harmonic  model  can  be  define  [11] 


X,  =A0  + £ a, 


COS 


2n  K(i  - 1)  n . 211  K(i  - 1) 

- — — ' + Pk  sin v — 


k=l  ' 


n 


n 


The  formula  for  determining  estimates  can  be  rewritten 
as  follows  using  Semi-amplitude  and  phase  angle  like 
[12-14]. 


rri 

X,  =A0  + £ v. 


sin 


k= 1 


2nK(i-l) 


+ < 


Where  Semi-amplitude  and  phase  angle  like  [12-14] 

V„=ilAl+Bl 

<pk  =tan"‘  (Ak/Bk)  k = l,2,3,....,m 

Tue  following  figures  give  as  example  of  sine  function 
which  is  a part  of  Fourier  model  the  main  aim  of  the 
research  was  to  build  w sine  ware  model  (equation)  using 
the  method  is  part  (1).  Then  use  the  build  model  to  graph 
any  sine  graph  [12-14]. 


Fig.l : example  for  sine  wave  function 


1MM| 

11^ 
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Let  us  have  the  following  control  point  is  periodic  in  lag  equal  to  22  and  we  want  to  estimate  the  Fourier  model  with  min 
error  [12-14]. 


Table.  1:  The  sine  wave  control  point 


1.1 

1.9 

2.7 

3.2 

4.0 

4.2 

5.0 

5.6 

5.9 

6.3 

6.3 

6.3 

6.2 

6.0 

5.8 

5.4 

5.0 

4.4 

4.2 

3.9 

3.2 

2.7 

2.2 

1.6  0.82 

0.25 

-0.4 

-1.01 

-1.5 

-1.9 

-2.25 

-2.6 

-2.6 

-2.7 

-2.5 

-1.9 

-1.65 

-1.20  -0.75 
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Then  by  using  Mini-Tab  package  we  write  a program  to  estimate  the  model  (the  Fourier  harmonic  model)  with  harmonic 
equal  to  4 and  by  using  the  method  in  part  1 we  estimate  the  model  with  following  parameter  [12-14]. 


K 

1 

2 

3 

4 


Semi-Amplitude 

4.38 

0.46 

0.12 

0.18 


Phase  angle 
-0.34 
0.11 
0.85 
0.23 


Where 


m n 

Ao=Z2>„/ 


i=i  j=i 


nxr 


Then  we  estimate  the  harmonic  function  from  the  observation  give  in  table  (1)  and  plot  it  with  actual  control  point  data  like 
in  figure  (2  ) the  result  show  the  model  was  good  representing  for  the  data  and  can  be  use  to  make  a sine  wave  plot  with 
small  error  [12-14]. 


Table.  2:  The  estimate  data  from  sine  model 


0.75 

1.77 

2.68 

3.41 

3.95 

4.35 

4.72 

5.11 

5.52 

5.96 

6.30 

6.41 

6.34 

6.14 

5.80 

5.43 

5.02 

4.70 

4.31 

3.82 

3.28 

2.75 

2.15 

1.61 

0.88 

0.25 

-0.40 

-1.1 

-1.61 

-1.92 

-2.41 

-2.46 

-2.51 

-2.54 

-2.43 

-2.01 

-1.77  -1.12  -0.254 


♦ Seriesl 


IV.  CONCLUSIONS 

From  the  research  we  conclude  the  following 

1-  The  sine  wave  function  can  be  built  using  a 
nonlinear  equation  or  a harmonic  function  with 
min  error  the  efficiency  of  the  model  increases 
if  the  number  of  observation  harmonic  number 
and  replication  number  Increase  but  this  will 
lead  to  increase  the  computational  time  in  the 
computer. 

2-  The  reconstruction  of  the  models  that  are 
basically  build  around  Sine  waves  makes  it 
possible  to  regenerate  the  original  forms  with 
new  features  added,  so  that  it  could  prevent 
further  counterfeiting. 


REFERENCES 

[1]  H.  Boche  and  B.  Farrell,  “Strong  divergence  of 
reconstruction  procedures  for  the  paley- wiener  space 
PW  7i  and  the  hardy  spaceHl,”  Journal  of 
Approximation  Theory,  vol.  183,  pp.  98-  117,  2014. 

[2]  J.  Higgins,  Sampling  theory  in  Fourier  and  signal 
analysis:  foundations.  Oxford  University  Press, 
1996. 

[3]  K.  Gr'ochenig,  Foundations  of  Time-Frequency 
Analysis.  Birkh  'auser,  2001. 

[4]  I.  Daubechies,  Ten  lectures  on  wavelets.  SIAM, 
1992,  vol.  61. 

[5]  C.  Shannon,  “Communication  in  the  presence  of 
noise,”  Proc.  IRE,  vol.  37,  no.  1,  pp.  10-21,  1949. 

[6]  D.  Gabor,  “Theory  of  communication,  part  1:  The 
analysis  of  information,”  Journal  of  the  Institution  of 
Electrical  Engineers -Part  III:  Radio  and 


www.ijaems.com 


Page  I 679 


International  Journal  of  Advanced  Engineering,  Management  and  Science  (IJAEMS) 
h ttps://dx.  doi.  org/1 0. 24001/iiaems.  3.6.10 

Communication  Engineering,  vol.  93,  no.  26,  pp. 

429-441,  1946. 

[7]  S.  Mallat,  “A  theory  for  multiresolution  signal 
decomposition:  the  wavelet  representation,”  Pattern 
Analysis  and  Machine  Intelligence,  IEEE 
Transactions  on,  vol.  11,  no.  7,  pp.  674-693,  1989. 

[8]  Generation  of  non-linear  curves  using  computer 
design  phd  thesis,  Firas  Husham  Al-Mukhtar 

[9]  Barnett. V.D  (1980)  fitting  straight  lines  the  linear 
functional  relationship  with  replicated  observations. 

[10] Sprent,  P.  (1981)  parallelism  and  occurrence  in 
linear  regression. 

[11] Enslein,K.,  Ralston  , A.  , & wilf  (1977)  statistical 
methods  for  digital  computers  , New  York  wiley. 

[12]  Draper  N.  R , & Smith  H.  (1967)  Applied 
Reqression  analysis  New  York  wiley 

[13]  Otens  R.  k and  Enochson  , L.D  (1972)  , Digital  time 
Series  Analysis  , New  York  , John  wiley  & Sons. 

[14] Otnes  , R.k  and  Enochson  L.D  (1978)  , Applied 
Time  Series  Analysis  , New  York  , John  wiley  & 

Sons. 


[Vol-3,  Issue-6,  Jun-  2017] 
ISSN:  2454-1311 


www.ijaems.com 


Page  I 680 


